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Abstract. The methodological approaches to compilation of landscape- 
geochemical maps for urban areas have been developed around cartograph- 
ic database which systematizes digital maps presenting functional zoning of 
the study area, its landscape structure and technogenic anomalies of heavy 
metals. The methods have been tested during GIS landscape-geochemical 
mapping of Eastern Administrative Okrug of the city of Moscow. A land- 
scape-geochemical map displaying the differentiation of urban landscapes 
in terms of conditions for heavy metal accumulation and according to con- 
tamination levels of soil and snow cover has been designed. 
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1. Introduction 

Landscape-geochemical mapping is one of the most important components 
of ecological mapping of urban areas. Its theoretical and methodological 
basis is formed by geochemistry of landscapes and envi ronmental geochem- 
istry (Saet et al. 1990, Perel'man & Kasimov 1999). Rapid development of 
GIS methods and technologies opens up new possibilities for compilation of 
landscape-geochemical maps and allows to benefit from systematic ap- 
proach in the analysis of various information on urban landscape structure 
and landscape functioning (Zhukovet al. 1999, Makarovetal. 2002, Kurba- 
tovaetal. 2004). 

In the course of research large amounts of data on components of urban 
landscapes are accumulated. These data require structuring and systemati- 
cal organization in order to improve efficient usability of the information. 
In addition there is a demand in providing information in the form of car- 
tographic coverage, because maps provide the most complete spatial char- 
acterization of the study objects. Such cartographic support can be present- 



ed by the cartographic database (CDB) containing a set of interrelated spa- 
tial data organized on certain rules and similar principles. 

The objectives of the study presented herein included: 

1) GIS mapping of megalopolis landscapes bysystematization of the data on 
landscape and functional structure in CDB. 

2) on the basis of the designed CDB, developing and testing the methods of 
compilation of synthetic landscape-geochemical map displaying the differ- 
entiation of urban landscapes according to the conditions of higher-priority 
pollutant behavior and defining environmental hazard of landscape pollu- 
tion in various functional zones of the city. 

As the object of the study a territory of Eastern Administrative Okrug (EAO) 
of the city of Moscow with a number of large industrial zones and highways, 
local heat station, combustion plant and some other pollution sources has 
been chosen. Among the major pollutants of urban landscapes heavy metals 
(HMs) come to the forefront because many of them are considered to be 
highlytoxic. 

The research on HMs in soils, snow and vegetation on the territory of EAO 
has been conducted since 1989, enabling to evaluate long-term dynamics 
and spatial structure of pollution in landscape components (Kasimov et al. 
2012, Labutina & Khaybrakhmanov 2012, Nikiforova et al. 2009-2011, 
Khaybrakhmanov 2011, Kasimov 1995). The results of the research provid- 
ed the basi s for the present study. 



2. The structure of CDB for landscape-geochemical 
studies 

Di gi tal maps form the content of CDB, therefore it's i important to determi ne 
the structure of attribute tables based on the types of both spatial (vector, 
vector-topological, raster) and graphic characteristics of objects (Luriye 
2008). As a preliminary stage it is necessary to develop the structure and 
content of the map support system def i ni ng the I ogi cal structure of CDB. 

The map support system was created for the territory of M oscow EAO. M ul- 
tispectral QuickBird satellite images (2009) with spatial resolution of 2.4 
meters, Ecological Atlas of Moscow (2000), cartographic and statistical 
Internet sources, data derived from geochemical soil and snow sampling 
served as the mai n sources of mappi ng. 

The created cartographi c materi al consi sts of the fol I owi ng bl ocks: 
1) Basic geographic information (base maps, satellite images); 



2) Maps of functional zones; 

3) M aps of the I andscape structure; 

4) Geochemical maps. 

Digital presentation of maps correspond to the established principles of GIS 
mapping including preliminary procedures for evaluating the quality of 
composed products, choosing symbolic systems, generalization, etc. 

The rules for establishing connections between individual attribute tables 
are determined by database management system (DBMS). I n our case, we 
used the DBMS ArcGIS 9.3 that provides management possibilities to deal 
with cartographic information. It defines the relational data model, which 
allows storing multiparameter spatial information in the form of related 
tables. Relations are defined by the columns of spatial objects' I D numbers 
on the map, which should be uniform for all interrelated objects. This ap- 
proach causes the importance of using the system of basic spatial cells for 
the storage of information. These eel Is were functional zones allocated with- 
in current research project when designing map support system. Cell num- 
bers are the same for different groups of indicators, spatially organized by 
functional zones. Eventually all the blocks of information were connected 
and a complete logical structure of CDB was developed (Fig. 1). This struc- 
ture connects the database blocks in the convenient logical model that al- 
lows joint consideration of the information to analyze and compile new 
maps, especi al ly compl ex and syntheti c ones. 




Figure L The logic structure of CDB 



3. The content of the database blocks 
3.1. Maps of functional zones 

The type of urban land use plays a leading role in the formation of techno- 
genic geochemical anomalies of pollutants (Methodological guidelines... 
2006), therefore in our research we developed the methods of functional 
zone mapping and compiled the map for EAO (Labutina & Khaybrakh- 
manov 2012, Khaybrakhmanov 2011). For the long time inventory data, 
various statistical data, topographical plans of cities have been the main 
sources of information. Updating the information, keeping it adequate was 
expensive and time-consuming. For the same reasons, aerial photography 
materials were not often used. Modern satellite images with spatial resolu- 
tion better than 5 m, being received with a several-day interval (GeoEye-1, 
WorldView-1, 1 konos, etc.) or even daily (RapidEye), have created the con- 
ditions for the development of remote sensi ng- based methods for the map- 
pi ng of urban areas (Aubrecht et al. 2009, Miller & Small 2003, Wilson et 
al.2003). 
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Figure 2. The map of functional zones for the territory of Moscow EAO 

Many urban objects (individual buildings, streets, driveways, railroad 
tracks, etc.) are confidently recognized on very high resolution satellite im- 
ages. Features and functions of objects are defined by a combination of sev- 
eral indicators - shape, size, color, relative position of objects, and also by 
logical deduction. At a lower resolution multispectral images indicate the 
properties of objects that are not visi ble on the black-and-white i mages. 



The map of functional zones for the territory of Moscow EAO has the scale 
150 000 (Fig. 2). The main data sources for this map are satellite images 
QuickBird (Digital Globe, USA) with a 0.61 m resolution in the panchro- 
matic band and 2.44 m - in multispectral band. The methods of automated 
and visual interpretation were used in the process of mapping. Green areas 
and water bodies were identified on the basis of computer image classifica- 
tion (Guo 2007), the structure of residential areas and city functional zones 
- by visual interpretation. 

The map of functional zones serves as a basis for revealing the factors of 
pollution and for the pollution assessment in the study area. During land- 
scape-geochemical studies such map helps to perform task-oriented sam- 
pling, to identify the patterns of wind-induced dust transfer according to 
building locations and building clustering in residential areas. Besides it 
makes possible evaluating the impact of pollution sources on landscapes. 

3.2. The maps of landscape structure 

During landscape-geochemical mapping of urban areas not only functional 
zoning but also displaying natural conditions and their anthropogenic 
transformation are important since these factors eventually define the be- 
havi or and fate of toxi cants and the dynami cs of I andscape pol I uti on . 

The major information on landscape structure of EAO was derived through 
the analysis and assessment of geological, geomorphologic, soil, vegetation, 
landscape maps from Ecological Atlas of Moscow (2000). In addition the 
maps designed by the authors were used (the maps of deposits, elementary 
geochemical landscapes, hydrous migration classes). All these maps were 
presented in a digital form in a separate block of the CDB devoted to land- 
scape structure and were considered together with the application of spatial 
analysis - overlay, morphometric operations, Gl S statistical analysis etc. As 
a result a set of maps characterizing landscape-geochemical relationships in 
urban landscapes was compiled. 

The map of deposits. Quaternary deposits in the study area are overlain 
by technogenic and cultural layers. Since such layers impact chemical com- 
position of soil cover and define initial levels of HM content they may con- 
tribute significantly to environmental and geochemical status of an urban 
landscape. For the analysis of deposit distribution over the territory of EAO 
the map displaying the differentiation of various (natural and technogenic) 
deposits in terms of their genesis, granulometry and thickness was de- 
signed. 
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Figure 3. The map of deposits i n EAO, the city of M oscow 

The map of elementary geochemical landscapes reveals catenary 
geochemical structure of landscapes and characterizes direction of lateral 
fluxes of polluting substances in soil cover of the urban area between inter- 
fluve positions and related depressions (Fig.4). As a basis of the map con- 
tent geochemical taxonomy of natural landscapes proposed by 
M.A.GIazovskaya was used (Glazovskaya 2002). Geochemical catenae be- 
tween elementary landscapes reflect the variety of soil and slope processes 
and form certain paragenetic groups characterized by certain types of mat- 
ter, energy and information exchange between subsystems. The map is a 
product of integrated and interrelated analysis of various maps - geomor- 



phologic, soil, landscape, a map of soil waterlogging and others. Taking into 
account the boundaries of local river watersheds 5 gradations (genera) of 
elementary landscapes have been distinguished on the map of EAO: eluvial, 
transdural, transaccumulative, superaquatic and aquatic. Due to levelled 
relief, small altitude differences and high anthropogenic transformation of 
relief features transeluvial landscapes on the territory of EAO are rather 
rare and eluvial landscapes often border superaquatic ones (Fig. 4). 




Figure 4. The map of elementary landscapes of EAO, the city of M oscow 

The map of hydrous migration of chemical elements. The intensity 
of hydrous migration of pollutants in soils and subsoils is controlled by re- 
dox and acid- base conditions. The combination of various redox and pH 
settings produces 12 hydrous migration classes (Perel'man & Kasimov 
1999). Anthropogenic impact on geochemical settings in soils of EAO, relat- 
ed to waterlogging processes and alkalinization, is shown on the map of 
hydrous migration of chemical elements. The map was compiled using the 
data collected during sampling in 2010 (Fig.5). In total the combination of 
redox and pH settings produced 7 hydrous migration classes which are 
shown on the map. 
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Figure 5. The map of the classes of hydrous migration of chemical elements in 
soi Is of EAO, the city of M oscow 

3.3. The maps of contamination of soils and snow cover with 
heavy metals 

Geochemical maps of this type display landscape-geochemical structure of 
the area under study. The levels of technogenic impact on the landscapes of 
EAO of the city of Moscow are characterized by the map of integral solid 
HM fallout onto snow cover (index Zd) and the map of integral soil pollu- 
tion Zc. For compilation of these maps the results of chemical analysis of 
snow and soil samples, collected in 2010 according to standard procedure 
(Methodological guidelines... 2006), were used. Geochemical sampling of 



snow and soil cover with the applied spacing (500-800m) made possible a 
quantitative assessment of recent levels of urban landscape contamination, 
identification of spatial trends in HM distributions and confirmation of 
technogeni c anomal i es contours. 

Soil samples and filtrates, containing sol id phase of snow, were analyzed for 
concentrations of 20 HMs by mass spectrometry and inductively coupled 
mass spectrometry (I CP- MS) at All-Russian Institute of Mineral RawMate- 
rials (VI M S) usi ng Perki n-EI mer devi ce, the USA. 

Dust load on landscapes (Pn ) was defined by concentrations of suspended 
matter in snow. The integral index of atmospheric addition (immission) of 
HMs was calculated using the formula Zd=£Kd - (n- 1), where Kd=DX/Df; 
DX, Df - fallouts of elements (mg/km2 per day) in urban and background 
environments respectively, n - the number of accumulating metals with Kd 
>15 (Methodological guidelines... 2006). 

Fallouts of certain elements onto snow surface were calculated using the 
formula D =Pn-C, where C is the concentration of an element in suspended 
matter (ppm). The levels of soil contamination with HMs were evaluated 
using integral index of pollution Zc = £Kc - (n- 1), where Kc=CX/Cf, CX, Cf 
- average concentrations of HMs (ppm) in urban and background soils re- 
spectively, n - the number of accumulating metals with Kc >1 (Saet et al. 
1990). 

The obtained results for each of 52 sampling sites were included into the 
CDB for spline interpolation (Geostatistical Analyst, ArcGIS 9.3) which en- 
abled to evaluate spatial patterns of integral HM pollution of snow cover 
and soils on the territory of EAO. As a result digital models of technogenic 
anomalies (Kasimov et al. 2012, Nikiforova et al. 2011) based on Zd and Zc 
indices were created. Their cartographic representation is given on Fig.6 
and 7. The analysis of these models resulted in identification of technogenic 
sources of HMs and isolation of HM associations on the basis of metal be- 
havior or accumulation intensity in the studied landscapes. 
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Figure 6. Integral pollution of snow cover with HMsin E AO of the city of Moscow 
according to addition (immission) index Zd 
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Figure7. Integral pollution of soils with HMsin EAOof the city of Moscow accord- 
ing to index Zc 



4. The application of CDB for compilation of synthetic 
landscape-geochemical maps for urban areas 

The database, comprising the blocks of interrelated digital maps, is de- 
signed to support integrated landscape-geochemical studies in urban lands. 
Such investigations suggest the compilation of synthetic landscape- 
geochemical maps which evaluate the differentiation of urban landscapes in 



terms of conditions of migration and accumulation of higher- priority pollu- 
tants. They also help to reveal environmental hazard of landscape pollution 
in various functional zones of the city. 

During compilation of large-scale landscape-geochemical maps for urban 
areas landscape-functional complexes as a major mapping units are used. 
The complexes are derived from nature and geochemical maps and from 
the map of functional zoning. Si nee functional zones are elementary spatial 
cells in theCDB, this can bedonethrough combination of digital map con- 
tours using the overlay procedure and automatic generalization in GIS 
software. 

Obviously in recent urban environment technogenic factors when compared 
to natural ones come to the forefront in terms of their impact on ecological 
situation (Zhukov et al. 1999, Perel'man & Kasimov 1999). It is confirmed 
by the values of rank correlation coefficients (r) estimated for the maps of 
soil and snow pollution compiled on the basis of Zc and Zd indices on one 
hand and the maps of functional zones and natural conditions on the other 
hand. Therefore during generalization of contours obtained through the 
overlay of the digital database maps, the contours were adjusted in accord- 
ance with those of functional zones. 

Consequently the polygons of 20 landscape- functional complexes which 
attri bute tables contain the information on landscape and functional struc- 
ture of the territory were obtai ned ( F i g.8) 

In the legend geochemical classification of urban landscapes is used as a 
basis (Perel'man & Kasimov 1999) and two major factors are taken into ac- 
count: 1- the intensity and the composition of technogenic loads defining 
the impact of major pollution sources, and 2 - landscape-geochemical set- 
tings, controlling accumulation and dispersion of the polluting substances 
and therefore determi ni ng the i impact effects. 

The legend is built in a matrix form. In the right-hand part of the legend the 
functional zones which define the addition of H Ms into landscapes are dis- 
played, and in the left-hand part landscape and natural features governing 
hydrous migration together with catenary patterns in HM distributions and 
the ability of HMs to accumulate at geochemical barriers are listed. The 
combination of these two parts enables to distinguish landscape- functional 
zones. Their features determine the formation and magnitude of HM geo- 
chemical anomalies in soils and snow (Fig.8). The shapes and the bounda- 
ries of the anomalies in different media are not uniform since chemical fea- 
tures of snow cover reflects the recent state of pollution whereas soil cover 
in urban environment accumulates pollutants during longer, many- years 
period. 
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Figure 8. Landscape-geochemical map of EAO of the city of M oscow 

On the designed map landscape-functional complexes are displayed by dif- 
ferent background colours. A colour itself characterizes the assignment of 
the functional zone, whereas colour intensity identifies the position of the 
elementary landscape within the catena: the more intensive colour is the 
more subordinate position the landscape occupies. Each taxon on the map 
is designated by an index which includes the name of the functional zone, 
the genus of the elementary landscape, the composition and genesis of de- 
posits and the thickness of technogenic layers. For example the index n- 
9i/i characterizes an eluvial landscape, located in an industrial zone, 



formed on fluvioglacial stony sands, overlain by cultural layer less than lm 
in thickness. 

The intensity of pollution of snow and soil covers with H Ms is displayed on 
the map by contour lines using integral indices of HMs immission Zd and 
soil pollution Zc, represented in geochemical block of maps in the CDB. To 
show pollution of snow vertical blue hatches are used. For soils red hatches 
with different spacing depending on the level of pollution are applied. The 
size and the colour of the hatched areas indicate cross location of techno- 
geni c anomal i es i n snow cover and soi Is. 



5. Conclusion 

On the basis of functional zoning and the analysis of landscape and geo- 
chemical structure of E AO in the city of Moscow the methods of specialized 
large-scale landscape-geochemical map compilation have been developed. 
Such maps display natural and anthropogenic factors and pollution levelsin 
depositing media. They show locations and sizes of technogenic geochemi- 
cal anomalies in urban landscape in relation to impact sources, landscape- 
geochemical settings and direction of lateral fluxes of polluting substances. 

The compilation of landscape-geochemical maps for urban areas is per- 
formed by assessment and generalization of considerable and various car- 
tographic materials which are arranged into interrelated blocks of digital 
maps in CDB. The database includes the map of functional zoning which is 
compiled on thebasisof high resolution satellite images and which contains 
a number of important characteristics necessary for environmental and 
geochemical investigations. Integrated GIS analysis of various maps of 
landscape structure enabled to reveal environmental features responsible 
for indication and deposition of HMs and other pollutants in urban land- 
scapes. Multi -elemental geochemical maps display spatial structure of 
tecnnogeni c anomal i es of H M s i n soi Is and snow cover of the studi ed area. 

The proposed methods have been tested with the example of urbanized eco- 
system located in south taiga - EAO of the city of M oscow. Using Gl S anal- 
ysis and generalization of various nature and nature- anthropogenic maps 
the landscape-geochemical map and its content have been designed. Its 
attribute table and the legend comprise natural and technogenic factors of 
HM accumulation-dispersion, levels of pollution and grades of ecological 
hazard i n terms of urban landscape pollution magnitude. 
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